Light-assisted ion-neutral reactive processes in the cold regime: 
radiative molecule formation vs. charge exchange 
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We present a combined experimental and theoretical study of cold reactive collisions between laser- 
cooled Ca + ions and Rb atoms in an ion-atom hybrid trap. We observe rich chemical dynamics which 
are interpreted in terms of non-adiabatic and radiative charge exchange as well as radiative molecule 
formation using high-level electronic structure calculations. We study the role of light-assisted processes 
and show that the efficiency of the dominant chemical pathways is considerably enhanced in excited 
reaction channels. Our results illustrate the importance of radiative and non-radiative processes for the 
cold chemistry occurring in ion-atom hybrid traps. 



Over the past few years, impressive progress has been 
achieved in the study of reactive collisions at ultralow en- 
ergies E co \\/kn <C 1 K HJ[2|. Recent landmark studies 
using neutral molecules highlighted the distinct quantum 
character of reactive processes in this regime and demon- 
strated new approaches for an unprecedented control of 
molecular collisions BHD. Ion-neutral reactions are an- 
other class of processes which exhibit different long-range 
interactions and therefore a different chemical behavior in 
comparison to neutrals II514TT11. With the development of 
hybrid traps in which laser-cooled atomic ions stored in a 
radiofrequency ion trap are combined with ultracold neu- 
tral atoms in a magneto-optical trap |[T2l - [T4l or a Bose- 
Einstein-condensate llT5l[T6l , the study of ion-neutral reac- 
tions in the energy range between 1 and 10 -3 Kelvin (usu- 
ally termed the "cold" regime) has recently become possi- 
ble. Under these conditions, only a few partial waves con- 
tribute to the collision so that resonance as well as radiative 
effects can become important ll5ll6l [T0l[T7ll . 

One key question pertains to the types of chemical pro- 
cesses which can occur in hybrid traps. So far, either fast 
near-resonant homonuclear charge exchange (in Yb-Yb + 
|[T3l0 or a slow loss of the atomic ions from the trap were 
observed (in Rb-Yb+ d and Rb-Ba+ HSI V For Rb-Yb+, 
the latter observation was rationalized in terms of radiative 
and non-radiative charge exchange [18]. The feasibility of 
molecular-ion formation has also been considered, and evi- 
dence for a radiative mechanism has recently been found in 
the Ca-Yb + system lfT4l . However, a general understand- 
ing of the interplay between these reactive processes and in 
particular the role of light remains to be established. 

In the current study, we present a combined experi- 
mental and theoretical study of ion-neutral reactive colli- 
sions in a Rb-Ca + hybrid trap. Our experimental results 
are interpreted using high-level electronic structure calcu- 
lations of the CaRb + potential energy curves (PECs) up 
to the twenty-second dissociation limit. We observe rich 
chemical dynamics which we rationalize in terms of non- 
adiabatic and radiative effects. We show that the efficiency 
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FIG. 1. (a) Schematic of the experimental setup, (b) Laser- 
cooling schemes for 40 Ca + and 8T Rb. (c) Superposed fluores- 
cence images of a Ca + ion Coulomb crystal (blue) and a cloud of 
ultracold Rb atoms (red) in the hybrid trap. 



of the dominant chemical processes (radiative molecule 
formation, radiative and non-radiative charge exchange) is 
considerably enhanced in excited reaction channels pop- 
ulated in the presence of radiation. Using Rb-Ca + as a 
model system, our results illustrate the reactive processes 
which can occur under the "cold" conditions of ion-atom 
hybrid traps. 

For the present study, an ion-atom hybrid trap was imple- 
mented by superimposing a linear radiofrequency ion trap 
for laser cooling 40 Ca + ions with a magneto-optical 
trap (MOT) fl]0 for 87 Rb atoms (Figure [l] (a)). Doppler 
laser cooling of the Ca + ions was achieved using two 
diode laser beams at 397 nm and 866 nm pumping on the 
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FIG. 2. (a) Series of Ca + ion laser-cooling fluorescence images 
as a function of the time of reaction with ultracold Rb atoms, (b) 
Pseudo-first-order analysis of the reaction kinetics of the experi- 
ment shown in (a). The uncertainties are smaller than the size of 
the symbols in the plot. 



transitions, respectively (Figure [I] (b)). Upon laser cooling, 
the ions form Coulomb crystals JH. The average ion ki- 
netic energies were dominated by the micromotion (the fast 
motion driven by the radiofrequency trapping fields) whose 
contribution was characterized by a comparison of experi- 
mental Coulomb crystal images with molecular-dynamics 
simulations (HO. In our experiments, the average colli- 
sion energies (E co u) were entirely governed by the ion ki- 
netic energies which were varied by changing the size and 
shape of the Coulomb crystals as discussed in Ref. 191 . 

The MOT was set up with two water-cooled solenoids in- 
stalled in vacuum to generate a quadrupolar magnetic field 
with a gradient of 20 G/cm. Laser beams around 780 nm 
in an optical-molasses configuration were used for cooling 
and repumping Rb on the (5s) 2 S 1 /2 — > (5p) 2 P^/2 transi- 
tion (see Fig. [T](b)). The number density and temperature 
of the Rb atoms in the MOT were established using stan- 
dard fluorescence measurement and time-of-flight meth- 
ods, respectively GUI . The fluorescence of either the ions 
or the neutral atoms was isolated using narrow-bandpass 
color filters and imaged onto a CCD camera (Figure [T](c)). 
A detailed description of the experimental apparatus and 
procedures will be given in a subsequent publication [21 j. 

In our experiment, the Ca + and Rb cooling lasers were 
alternately blocked using a mechanical chopper at a fre- 
quency of 1000 Hz in order to prevent photoionization of 
Rb out of the (5p) 2 P^/2 level by 397 nm photons. Un- 
der these conditions, the typical Rb number densities and 
temperatures achieved were on the order of 1 x 10 9 cm -3 
and T — 150 — 200 fiK, respectively. The chopping of the 
cooling-laser beams only negligibly affected the kinetic en- 
ergies of the ions and their overlap with the Rb cloud. 

Reactive collisions between Ca + and Rb lead to a de- 
crease in the number of Ca + ions in the Coulomb crys- 
tals as shown in Figure [2] (a). When both species are 
laser cooled, the populations are distributed over the 
(4s) 2 5i /2j (4p) 2 Pi /2 and (3d) 2 D 3/2 states of Ca+ and 
the (5s) 2 Si/2 and (bp) 2 P$/2 states of Rb (compare Fig. 
[T] (b)). Reactive collisions occur in excited states of the 
Rb-Ca + system and the observed reaction rates represent 
an average over all possible channels. Note that simulta- 
neous excitation of both species does not occur because of 
the alternate chopping of the cooling laser beams. 

Rate coefficients were determined by measuring the de- 
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FIG. 3. (a) Variation of reaction rate constants k with Ca + level 
populations. The bars indicate the relative level populations ob- 
tained by varying the detunings of the cooling-laser beams, their 
intercepts with the y axis give the values of the corresponding 
rate constants. Each data point corresponds to an average of 
three consecutive measurements (average statistical uncertainty 
(2a) Ak = 3 x 10" 12 cm 3 s" 1 ). (b) Rate constant as a func- 
tion of the average collision energy {E co \\)/k^. The error bars 
denote the statistical uncertainty (2a) of the measurements, (c) 
Resonant-excitation mass spectra of Coulomb crystals (i) before 
reaction, (ii) after reaction. The dashed vertical lines indicate 
single-ion motional frequencies. 



crease of the Ca + Coulomb crystal volume V as a func- 
tion of the reaction time and fitting the results to a pseudo- 
first-order rate expression \n(V/Vo) = k't Q where 
V denotes the initial crystal volume and t is the reac- 
tion time, see Figure [2] (b). Second-order rate coeffi- 
cients k were obtained by dividing the pseudo first-order 
rate constants k' by the average density of Rb atoms 7V Rb 
in the MOT: k — k' /N^. For the measurement dis- 
played in Figure[2j the rate coefficient was established to be 
k — 2.5(9) x 10 -11 cm 3 s _1 , only two orders of magni- 
tude smaller than the collisional (Langevin) rate coefficient 
(k L = 3 x 1(T 9 cm 3 s- 1 d). 

When the Ca + cooling lasers were switched off, a rate 
coefficient k s = 3(1) x 10 -12 cm 3 s _1 was obtained. This 
value proved to be insensitive to the excited-state popula- 
tion of Rb which was varied in the range from 2% to 5% 
compatible with a stable operation of the MOT. However, 
the ions are likely to heat up in the absence of laser cooling 
which also reduces their overlap with the cold-atom cloud. 
Therefore, an accurate value for collision energies cannot 
be given under these conditions and this value for k s must 
be regarded as an estimate of the rate coefficient in the low- 
est reaction channel Ca + (4s)+Rb(5s). 

To assess the contribution of the excited Ca + channels, 
we modified the populations in the Ca + states by vary- 
ing the frequency detuning of the 397 nm and 866 nm 
lasers in the range 20-100 MHz and 0-70 MHz, respec- 
tively. In these measurements, it was ensured that the ion 
cloud always remained Coulomb-crystallized so that the ef- 
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feet of the detuning on the average ion energies is negligi- 
ble. Figure[3](a) shows the resulting rate coefficients which 
clearly correlate with the population in the Ca + (4p) 2 Pi/ 2 
state, whereas no correlation with the populations in the 
(4s) 2 Si/2 and (3d) 2 D 3 / 2 levels is observed. The popu- 
lations were inferred from an Einstein rate-equation model 
of the laser excitations in the Ca + three-level system, with 
an estimated relative uncertainty of 20% [21 j. 

We fitted the data in Fig. [3] (a) to a kinetic model ac- 
counting for the relevant reaction channels. Taking into ac- 
count the chopping of the Ca + cooling laser beam, the mea- 
sured rate coefficient can be expressed as k = \ [k s (p s + 
1 ) + k p p p + k d p d ] . Here, k s , k p , k d and p s , p p , p d stand for 
the rate coefficients and populations in the Ca + (4s), (4p) 
and (3d) levels. We find k p = 1.5(6) x 1(T 10 cm 3 s" 1 
andfc a , d < 3 x KT 12 cm 3 s" 1 . 

Rate coefficients as a function of the collision en- 
ergy are displayed in Fig. [3] (b). Within the range 
(E co \x)/k# =200 mK-20 K, we observe an increase by a 
factor of two. A detailed analysis of these observations 
will be published in a subsequent article [21 j. 

The chemical identity of the reaction products was estab- 
lished using resonant-excitation mass spectrometry of the 
Coulomb crystals [23]. The spectra taken before immers- 
ing the Ca + Coulomb crystal in the MOT only show a sin- 
gle resonance corresponding to the excitation of the Ca + 
ions (Figure [3] (c) (i)). By contrast, the spectra recorded 
after the reaction (Figure |3](c) (ii)) exhibit in total four dis- 
tinct resonances which are identified as Ca + , Rb + , CaRb + 
and RbJ by comparison with the expected single-ion exci- 
tation frequencies (indicated by dashed vertical lines). 

Whereas the Rb + product is a result of charge exchange 
between Ca + and Rb, CaRb + can only be formed in an 
association reaction whereby the collision complex is sta- 
bilized either collisionally or by the spontaneous emis- 
sion of a photon (radiative association, RA). Collisional 
stabilization is negligible at the low Rb densities in the 
MOT. Moreover, RbJ is only observed in the presence of 
CaRb + indicating that this product results from the reaction 
CaRb + +Rb— ^ RbJ + Ca which is exothermic by ^0.4 eV 
according to our calculations. Product-ion peaks for Rb + 
and CaRb + were observed over the range of cooling laser 
detunings used in the present study and when the Ca + cool- 
ing laser was switched off, indicating that these ions are 
produced in the ground as well as in the excited reaction 
channels l2T1l . 

Figure [4] displays the RbCa + potential energy curves 
(without spin-orbit interaction) up to the twenty-second 
dissociation threshold Rb(5s)+Ca + (4p) which can be 
reached in the present light-induced dynamics. Currently, 
the only reliable way to describe such highly-excited states 
consists in treating explicitly only the two valence electrons 
of CaRb + in the field of two ionic cores Rb + and Ca 2+ rep- 
resented by effective core potentials J24l[25]]. Correlation 
between core and valence electrons is modeled via effective 
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FIG. 4. Computed RbCa + potential energy curves (PECs) in 
the regions of (a) Rb(5s)+Ca + (4s), (b) Rb(5s)+Ca + (3d), (c) 
Rb(5s)+Ca + (4p). The PECs of the relevant entrance channels 
of the reaction are highlighted with thick lines. PECs for 1 A and 
3 A symmetries are omitted for clarity. Downward arrows suggest 
possible pathways for formation of ground-state RbCa + ions by 
radiative association. In (c), only PECs of 1 S + symmetry are 
displayed for clarity. 



core polarization potentials. Therefore, a full configuration 
interaction (FCI) calculation is achievable in a configura- 
tion space built from a large basis set of Gaussian orbitals 
for the valence electrons (23 [26). The complete set of cal- 
culations will be presented elsewhere EH . 

Over the energy range of about 41000 cm -1 , two kinds 
of asymptotic channels are present: a Rb (5s or bp) atom 
colliding with a Ca + (4s, 3d, or 4p) ion (4 channels), and a 
Rb + ion colliding with a Ca atom (18 channels). In our ap- 
proach, the level energies of Rb and Ca + are constrained to 
the experimental values. The energies of the effective two- 
electron Ca atom are obtained through the FCI calculation, 
with an accuracy of a few hundred wave numbers typical 
for such calculations. Therefore, in the congested region 
around the Rb(5s)+Ca + (4p) asymptote (Figure [4] (c)) the 
order of several asymptotic limits is changed, preventing a 
detailed state-to-state modeling of the dynamics. Neverthe- 
less, the calculations provide a useful guide to a qualitative 
understanding of the underlying reaction mechanisms. 

Figure [4] (a) shows the region already studied theoreti- 
cally in Ref. (22]]. Note that the Rb(5s)+Ca + (4s) entrance 
channel is not the lowest dissociation limit. Non-radiative 
charge exchange (NRCE) induced by non-adiabatic cou- 
plings can occur around the crossings with the (1) 3 II state 
correlating with the Rb + +Ca(4s4p 3 P) asymptote [22]. 
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Moreover, radiative relaxation from the (2) 1 S + entrance 
channel to the state can either lead to radiative 

charge exchange (RCE) or RA forming molecular ions pro- 
vided that the relevant Franck-Condon (FC) factors are suf- 
ficiently large. The dashed arrow (i) in Figure [4] (a) sug- 
gests a transition for which favorable FC overlap for RA 
can be expected. This picture is compatible with our ex- 
perimental observation of both RA and CE products in this 
channel. We note that the upper limit for the experimental 
rate coefficient k s < 3 x 10 -12 cm 3 s _1 determined above 
is close to value of « 1 x 10 -12 cm 3 s _1 estimated from 
theoretical NRCE cross sections reported in Ref. [22 J . 

In the case of the Rb(5s)+Ca + (3d) entrance channel, 
no close-lying molecular state asymptotically correlating 
with neutral Ca exist which would enable efficient NRCE 
(Figure [4] (b)). This is in agreement with our experi- 
mental finding of a small reaction rate for this channel 
(k d < 3 x 1CT 12 cm 3 s" 1 ). RCE and RA connecting to 
lower-lying states are in principle possible in various sym- 
metries. However, inspection of the relevant PECs suggests 
that the FC factors are small. 

Figure [4] (c) shows the complicated network of the 
PECs for the 1 S + states in the region around the 
Rb(5s)+Ca + (4£>) entrance channel resulting from the 
large number of dissociation limits within a small energy 
range. The situation is comparable for the other sym- 
metries which are not shown for clarity. From Figure 
[4] (c), a mechanism explaining the large reaction rates 
in this channel becomes apparent. At large distances, 
the Rb(5s)+Ca + (4£>) PEC is attractive (varying as -R~ A ) 
and isolated from other channels with the same asymp- 
totic charge state, e.g., Rb(5s)+Ca + (4s, 3d). By contrast, 
while the long-range PECs of the states correlating with 
the Rb + +Ca* asymptotes are also attractive, they soon turn 
into repulsive curves as being Rydberg states converging 
to the Rb + +Ca + ionization threshold. As a consequence, 
the PEC of the entrance channel (the (13) 1 S + state) under- 
goes avoided crossings with lower curves correlating with 
charge-transfer asymptotes. Non-adiabatic transitions at 
these crossings will result in a large rate for NRCE lead- 
ing to the formation of Rb + . Moreover, through consec- 
utive non-adiabatic transitions in the intricate network of 
avoided crossings visible in Fig. [4] (c), the collision part- 
ners can approach to short internuclear distances where 
many channels are open for RCE and RA in all symme- 
tries. A possible radiative relaxation pathway to the abso- 
lute ground state is suggested by the downward arrow (ii). 
Additionally, radiative relaxation at long range can also 
contribute to the observed charge-exchange and molecule- 
formation processes. 

This mechanism for radiative association and charge 
transfer from excited states is reminiscent of the pro- 
cesses leading to "radiative escape" of ultracold atoms 
from MOTs IT7I : excitation occurs in the asymptotic re- 
gion in the early stages of the collision, followed by ra- 
diative relaxation at shorter range. In general, the multi- 



tude of different charge-transfer channels in ion-atom hy- 
brid systems leads to a considerable increase of the density 
of states in comparison to purely neutral collision systems. 
This situation implies a higher probability for favorable 
Franck-Condon overlaps which increases the possibilities 
(and therefore the rates) for RA and RCE out of excited 
channels. This situation can be generalized to all alkaline 
atom-alkaline earth ion systems and it can thus be expected 
that RA and RCE are generally important processes in the 
collisional dynamics of these species. 

In summary, our results illustrate the rich chemical 
dynamics exhibited by even the simplest possible (i.e., 
atomic) ion-neutral collision systems in the cold regime. 
In particular, we highlighted the complex interplay be- 
tween non-adiabatic effects, radiative charge exchange and 
molecule formation. Our findings underline the importance 
of light in enhancing these processes and suggest that effi- 
cient photoassociation of molecular ions using red detuned 
laser radiation might be feasible in the Rb-Ca + system. 
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